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Abstract. Gutzwiller's trace formula for the semiclassical density of states in 

a chaotic system diverges near bifurcations of periodic orbits, where it must be 

replaced with uniform approximations. It is well known that, when applying these 

f*^ ' approximations, complex predecessors of orbits created in the bifurcation ("ghost 

^^ ■ orbits" ) can produce pronounced signatures in the semiclassical spectra in the vicinity 

of the bifurcation. It is the purpose of this paper to demonstrate that these ghost orbits 

themselves can undergo bifurcations, resulting in complex, nongeneric bifurcation 

Q\ . scenarios. We do so by studying an example taken from the Diamagnetic Kepler 

^h ' Problem, viz. the period quadrupling of the balloon orbit. By application of normal 

^■y form theory we construct an analytic description of the complete bifurcation scenario, 

which is then used to calculate the pertinent uniform approximation. The ghost orbit 

Q ' bifurcation turns out to produce signatures in the semiclassical spectrum in much the 

same way as a bifurcation of real orbits would. 
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1. Introduction 

Since its discovery in the early 1970s, Gutzwiller's trace formula Jl|, H has become a 
widely used tool for the interpretation of quantum mechanical spectra of systems whose 
classical counterpart exhibits chaotic behaviour. It represents the density of states of 
the quantum system as a sum over a smooth part and fluctuations from all periodic 
orbits of the classical system, where the contribution of a single periodic orbit reads 

y giCSpo/Ti-fMpo) 

v /|det(M po -/)| 

with T po , S po , M po , fi po denoting the orbital period, action, monodromy matrix, and 
M aslov index, respectively. This formula assumes that all periodic orbits can be regarded 
as isolated, which is the case, in particular, for completely hyperbolic systems. In the 
generic case of mixed regular-chaotic dynamics, however, the formula fails whenever 
bifurcations of periodic orbits occur, because close to a bifurcation periodic orbits 
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approach one another arbitrarily closely. The failure of the formula manifests itself 
in divergences of the periodic orbit contributions (|1|). 

The generic cases of period-m-tupling bifurcations were studied by Ozorio de 
Almeida and Hannay || [| , who derived uniform semiclassical approximations by taking 
into account all orbits involved in a bifurcation collectively. Their solutions were refined 
by Sieber and Schomerus [EL ^ [7|, who derived uniform approximations for all types 
of bifurcations of codimension one in generic Hamiltonian systems with two degrees of 
freedom. Their formulas smooth the divergence in Gutzwiller's trace formula, and, in 
contrast to the approximations in ||, || , asymptotically approach the result of the trace 
formula ([!]) for isolated periodic orbit contributions as the distance from the bifurcation 
increases. As a characteristic feature, uniform approximations require the inclusion 
of complex "ghost orbits". At the bifurcation points, new periodic orbits are born. 
However, before they come into being, the orbits possess predecessors - ghost orbits 
- in the complexified phase space. As was shown by Kus et al. ||, some of these 
ghost orbits, which in the limit h — > yield exponentially small contributions, have to 
be included in Gutzwiller's trace formula ([!]). As a result, in constructing a uniform 
approximation complete information about the bifurcation scenario including the ghost 
orbits is required. 

A closer inspection of the bifurcation scenarios encountered in practical applications 
of uniform approximations reveals that bifurcations of codimension two, although they 
cannot generically be observed if only one control parameter is varied, can nevertheless 
have an effect on semiclassical spectra, because in their neighbourhood two bifurcations 
of codimension one come close to each other, and therefore have to be treated collectively. 
Examples of that situation have been studied by Main and Wunner |9|, [[(J as well as 
Schomerus and Haake [O, IT2 ]. 

The bifurcation scenarios described in the literature so far involve bifurcations 
of real orbits only. However, one should expect bifurcations of ghost orbits also to 
be possible and of particular importance for complicated bifurcation scenarios with 
codimension greater than one. It is the purpose of this paper to demonstrate that ghost 
orbit bifurcations do indeed occur and have a pronounced effect on semiclassical spectra. 
To this end, we present an example taken from the Diamagnetic Kepler Problem. It 
turns out that even the analysis of the period-quadrupling bifurcation of one of the 
shortest periodic orbits in that system requires the inclusion of a ghost bifurcation. 

The appearance of ghost orbit bifurcations represents an additional challenge for 
the construction of uniform approximations. It will turn out that normal form theory 
allows treating both ghost bifurcations and bifurcations of real orbits on an equal 
footing. Consequently, ghost bifurcations contribute to uniform approximations in much 
the same way as real bifurcations do. From these observations we conclude that the 
occurence of ghost bifurcations in systems with mixed regular-chaotic dynamics is not 
a very exotic but rather quite a common phenomenon. 

The organization of the paper is as follows: In section ^| we describe the bifurcation 
scenario of the example chosen in detail. Section ||] provides the general form of uniform 
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approximations. Section ^ presents the normal form describing the bifurcation scenario 
in point and the discussion of how the ghost bifurcation can be included in the normal 
form. In section [| we determine the normal form parameters to quantitatively describe 
the bifurcations, and in section |6] the uniform approximation is evaluated. 

2. The bifurcation scenario 

As an example, we study the hydrogen atom in a magnetic field, which has been 



described in detail, e.g., in Refs. [13], [14], [Lql . We assume the nucleus fixed and regard 



the electron as a structureless point charge. If the magnetic field is directed along the 
z-axis, the nonrelativistic Hamiltonian describing the electron motion reads [in atomic 
units, with 7 = 5/(2.35 x 10 5 T) the magnetic field strength] 

H= 1 -p>+ 1 -lL z+ 1 - 1 >( X i + y*)- 1 - = E. (2) 

Here, r is the distance from the nucleus, and L z denotes the angular momentum along 
the field axis, which is conserved because of the rotational symmetry around that axis. 
In the following we restrict ourselves to the case where L z = 0. 

To further simplify the Hamiltonian, we exploit its scaling property with respect to 
the magnetic field strength 7. In scaled coordinates and momenta 

r = 7 2 / 3 r, p = 7~ 1/3 p (3) 

the Hamiltonian assumes the form 

H = r 2/ "H = \tf + \^ + tf)-\ = E. (4) 

Thus, the classical dynamics does not depend on the energy E and field strength 7 
separately, but only on the scaled energy E = r )~ 2 / z E. From the scaling prescriptions 
(H) and (f|) we derive the scaling laws for classical actions and times as 

S = 7 1 / 3 S f = 1 T . (5) 

Due to the Coulomb potential, the Hamiltonian (Q) is singular at f = 0. The 
equations of motion can be regularized by introducing semiparabolical coordinates 

/i 2 = f + z , v 2 = r — z , (6) 

and a new time parameter r defined by alt = 1r dr. Finally, the regularized Hamiltonian 
is obtained as [O, [14], |15[ 



H 2 



\ (?l + Pt)-E (a 2 + V) + i„V [f + V) = 2 , (7) 

and Hamilton's equations of motion read [with primes denoting derivatives d/dr], 

t*' = Pv » Pv = 2 ^n - -/iz/ 2 (2/! 2 + v 2 ) , 

{ (8) 

v' = Pu , p' v = 2Eu - -fi 2 u(fi 2 + 2z/ 2 ) . 

These equations are free of singularities and can easily be integrated numerically. When 
using them, we must keep in mind that the definition (Q) determines the semiparabolical 
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coordinates /i, v up to a choice of sign only, giving a many-to-one coordinate system. 
Thus, if we integrate the equations of motion (|8|) until the trajectory closes in (fi, v)- 
coordinates this may correspond to more than one period in the original configuration 
space. Furthermore, we have to identify orbits which can be transformed into one 
another by reflections at the coordinate axes. 

We now complexify our phase space by allowing coordinates and momenta to 
assume complex values. This extension allows us to look for ghost predecessors of 
real orbits born in a bifurcation |§. 

At any given scaled energy E, there is a periodic orbit along the magnetic field 
axis. For sufficiently low negative E it is stable, while, as E /* 0, it loses and regains 



stability infinitely often [1G, [L3L Stability is lost, for the first time, at E = —0.391. At 



this scaled energy, the so-called balloon orbit |17j is born as a new, stable periodic orbit. 
As the scaled energy increases further, the orbit exhibits all kinds of period-m-tupling 
bifurcations before finally turning unstable at E = —0.291. 

Here, we consider the period-quadrupling bifurcation of the balloon orbit at E c = 
—0.342025. For E > E c , two real satellite orbits of quadruple period exist. As E \ E c , 
they collide with the balloon orbit and form an island-chain bifurcation as described in 
|3|, f4|. The real orbits are shown in figure [I] at scaled energy E = —0.340. The solid 
and dashed curves represent the stable and unstable satellite orbits, respectively. For 
comparison, the balloon orbit is shown as a dotted curve. Below E c , no real satellite 
orbits exist. Instead, there are a stable and an unstable complex ghost satellite. The real 
and imaginary parts of the stable and unstable ghost orbits at scaled energy E = —0.343 
are drawn as solid and dotted curves in figure ^. Note that the imaginary parts are small 
compared to the real parts. As predicted by normal form theory (see section |4.1j ), both 
satellites coincide with their complex conjugates, whence the total number of orbits is 
conserved in the bifurcation. 

The orbits described so far form a generic type of period-quadrupling bifurcation as 
investigated by Sieber and Schomerus f7j. However, the classical periodic orbit search in 
the complexified phase space reveals the existence of an additional ghost orbit at scaled 
energies around E c . The shape of this orbit is shown as the dashed curve in figure 
0. It is similar to the stable ghost satellite originating from the period-quadrupling 
of the balloon orbit. When following this ghost orbit to lower energies we find another 
bifurcation at E' c = —0.343605, i.e., slightly below the bifurcation point E c = —0.342025 
of the period-quadrupling of the balloon orbit. At energy E = E' c the additional ghost 
orbit (dashed curve in figure 0) collides with the stable ghost satellite of the period- 
quadrupling bifurcation (solid curve in figure H), and these two orbits turn into a pair of 
complex conjugate ghost orbits. Their shapes are presented at scaled energy E = —0.344 
by the solid and dashed curves in figure |3|. The imaginary parts clearly exhibit the loss of 
conjugation symmetry described, if the aforementioned symmetry of the semiparabolical 
coordinate system is taken into account. The dotted curve in figure |3] is the unstable 
ghost satellite which already exists at higher energy E > E' c (see the dotted curve in 
figure |). 
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It is important to note that the second bifurcation at E = E' c involves ghost 
orbits only. This kind of bifurcation has not been described in the literature so 



far; in particular, Meyer's classification of codimension-one bifurcations [IS] contains 
bifurcations of real orbits only and does not cover ghost bifurcations. The existence of 
the ghost orbit bifurcation implies that the results of Ref. |7| for the uniform semiclassical 
approximation for the generic period-quadrupling bifurcation cannot be applied to the 
more complicated bifurcation scenario considered here. As in cases described before by 
Main and Wunner || |l(| and Schomerus and Haake [|ll|, [l2j , the closeness of the two 



bifurcations requires the construction of a uniform approximation taking into account 
all orbits involved in the successive bifurcations collectively. Thus, the ghost bifurcation 
at E' c turns out to contribute to the semiclassical approximation in the same way as 
a real bifurcation would, as long as we do not go to the extreme semiclassical domain 
where the bifurcations can be regarded as isolated. We will demonstrate in section f| 
that the techniques of normal form theory can be extended as to include the description 
of ghost bifurcations. 

The construction of the uniform approximation requires the knowledge of the 
periodic orbit parameters of all orbits participating in the bifurcation scenario. The 
numerically calculated parameters are shown in figure f| as functions of the scaled energy 
E. Part (a) of this figure displays the actions of the periodic orbits, where the action 
of four repetitions of the balloon orbit was chosen as a reference level (AS = 0). This 
kind of presentation exhibits the sequence of bifurcations more clearly than a plot of the 
actual action integrals. Around E c , we recognize two almost parabolic curves, indicating 
the actions of the stable (upper curve) and unstable (lower curve) satellites, respectively. 
At E c , the curves change from solid to dashed as the satellites become complex. Below 
E c , the unstable ghost satellite does not undergo any further bifurcations in the energy 
range shown, whereas the stable satellite collides, at E' c , with the additional ghost orbit, 
which can clearly be seen not to be involved in the bifurcation at E c . Below E' c , these 
two orbits are complex conjugates of each other. Thus, the real parts of their actions 
coincide, whereas the imaginary parts are different from zero and have opposite signs. 

Analogously, figure f|b displays the orbital periods. Here, no differences were taken, 
so that the fourth repetition of the balloon orbit, which is always real, shows up as a 
nearly horizontal solid line at T «s 5.84. The other orbits can be identified with the 
help of the bifurcations they undergo, in the same way as discussed above. Finally, 
figure |]c presents the traces of monodromy matrices minus two. These quantities agree 
with det(M — /) for systems with two degrees of freedom. At E c and E' c , they can be 
seen to vanish for the orbits involved in the bifurcations, causing the divergences of the 
periodic orbit amplitudes ([[]) at the bifurcation points. 

The generic period-quadrupling bifurcation at E c can be described with the help 
of lowest order normal form theory presented in section (4.1[ As the ghost bifurcation 
is approached, this description fails. However, the influence of the additional ghost 
orbit can be taken into account by including higher order terms in the normal form as 



discussed in section 4.2. 
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3. The general form of the uniform approximation 

Before we return to classical normal form theory in section |], we introduce, in 
this section, the basic formulas for the quantum density of states necessary for the 
construction of the uniform semiclassical approximation. The density of states of a 
quantum system with the Hamiltonian H can be expressed with the help of the Green's 
function G{E) = (E - H)- 1 as 

d(E) = — Im Tr G(E) , (9) 

7T 

where the trace of the Green's function can be evaluated in the coordinate representa- 
tion, 

Tr G(E) = f dx'dx5(x' - x)G(x'x, E) . (10) 

The basic steps in the formulation of periodic orbit theory Jj], are to replace the Green's 
function G(x'x, E) with its semiclassical Van Vleck-Gutzwiller approximation and to 
carry out the integrals in stationary-phase approximation. For systems with two degrees 
of freedom the semiclassical approximation to the Green's function reads 

G&x,E) = —j== £ s/Dexp^S(x'x,E)-fy} . (11) 

Here, the sum extends over all classical trajectories running from x to x' at energy 
E, S is the action of a trajectory, v its Maslov index, and D is defined by the second 
derivatives of the action, 

/ d 2 s d 2 s \ 
£> = det <»f d fl). (12) 

\ dEdx dE 2 J 

The contribution of a single orbit to the density of states can be evaluated by introducing 
coordinates parallel and perpendicular to the orbit. The integration along the orbit can 
then be performed in a straightforward fashion. Finally, Gutzwiller's trace formula (]T|) 
for isolated periodic orbits is obtained by integrating over the coordinates perpendicular 
to the trajectory using the stationary-phase approximation. It is this last step which 
fails close to a bifurcation, where periodic orbits are not isolated. 

To find an expression for the density of states which is valid close to bifurcations, it is 
convenient to go over to a coordinate-momentum representation of the Green's function. 
Close to a period-m-tupling bifurcation of a real orbit the uniform approximation takes 
the form B 



d{E) = ^nl? Re I dydp ydE\\ dyW y exp {]i( § + yp 'y)- i r} ' (13) 

with y, p' y , and S defined as follows. Let (y, p y ) be the canonical variables in the Poincare 
surface of section perpendicular to the bifurcating orbit, and (y',p' y ) the corresponding 
variables after the period-m cycle of the orbit. The function 

S(p' y y,E) = S(y'y,E)-y'p' y (14) 
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is the Legendre transform of the action integral with respect to the final coordinate 
or, in other words, the generating function of the Poincare map for m periods of the 
bifurcating orbit in a coordinate-momentum representation. 
The function in the exponent in (p~3|), 

f(y,Py,E) = S(y,Py,E)+yPy, (15) 

is stationary at the fixed points of the m-traversal Poincare map, that is, fixed points 
of / correspond to classical periodic orbits. In the spirit of catastrophe theory, we now 
relate / to a given ansatz function $, which has the same distribution of stationary 
points, by a smooth invertible change of coordinates ip as 

f(y,p' y ,E) = S (E)+m}(y,p' y ;E);E) . (16) 

Here, ip is assumed to keep the origin fixed, xjj(0, 0; E) = (0, 0), and Sq(E) is the action 
of the central bifurcating orbit. Inserting this ansatz, we obtain 
1 



d ^=2^? ite «p{^- i H 




ex P |^$(r,p;; 



;i7) 



where Hess denotes the Hessian matrix and the coefficient 

| Hess $ | 
\ |Hess/| 

is unknown. At a stationary point of $, it can be shown to assume the value 

{m}T 



X 



sp 



TrM-2| 



Hess $ I 



(19) 



where T and M denote the period and the monodromy matrix of the corresponding 
classical orbit, and the notation {m} is meant to indicate that this factor does not 
occur at the satellite orbits. 



4. Normal-form description of the bifurcation 

To evaluate the uniform approximation (|17|), we need to find a suitable ansatz function 



$. Normal form theory as developed by Birkhoff |19| and Gustavson |2(| provides us with 
a systematic way to construct such ansatz functions. We will first investigate the lowest 
nontrivial order of the normal form expansion which describes generic bifurcations. 
Then, we will show that higher-order terms in the expansion can account for the more 
complicated bifurcation scenario studied here. 



4-1. The generic period- quadrupling bifurcation 

The simplest normal form describing the generic period-quadrupling bifurcations reads 



$ = £ I + al 2 + hi 2 cos(4<^) . 



(20) 
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This normal form is expressed in terms of canonical (action-angle) polar coordinates 
(I,(f), which are connected to Cartesian coordinates (p, q) by 

p=v2Jcosv?, g = v2/sin<y9. (21) 

To establish the connection with classical periodic orbits, we need to determine 
the stationary points of the normal form and then evaluate its stationary values. The 
stationary points are given by the equations 

= -^ = -4J 2 6sin(4^) , 

d(p (22) 

= — = e + 2al + 2bl cosUcp) . 
ol 

The first of these equations yields 

sin(4<z>) = , 

V ^ ! (23) 

cos(4</>) = a = ±1 , 

so that the second equation reads 
e + 2(o + ab)I = , 
with its solution 

(24) 



2(o + ab) 

In addition, there is the central periodic orbit at I = 0, which does not show up as a 
stationary point because the polar coordinate system ( pT|) is singular at the origin. 

To interprete these results, we observe that according to its definition (^1|) the 
coordinate I is positive for real orbits and that the action difference $(/, ip) is real 
for real /, (p. Therefore, negative real solutions I correspond to ghost orbits which are 
symmetric with respect to complex conjugation and thus have real actions, whereas a 
complex I indicates an asymmetric ghost orbit. 

Now, if \a\ > \b\, the two solutions I± have the same sign, which changes at e — 0. 
So, the two satellite orbits change from two real orbits to two ghosts at the bifurcation 
point e = 0, forming an island-chain-bifurcation. If \a\ < \b\, I + and /_ have different 
signs, so that on either side of the bifurcation there is one real and one ghost satellite. At 
e = 0, the satellites collide with the central orbit, forming a touch-and-go-bifurcation. 

4-2. Generalization to nongeneric bifurcations 

To describe a sequence of two bifurcations, we need to include higher-order terms in the 
normal form. Here, we adopt the normal form 

$ = el + al 2 + bl 2 cos(4</?) + c/ 3 (l + cos(4y?)) (25) 

given by Schomerus PTJ as a variant of the normal form used by Sadovskii and Delos 



to describe a sequence of bifurcations close to a period-quadrupling. It turns out to 
qualitatively describe the sequence of bifurcations encountered here for suitably chosen 
parameter values. 
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The stationary points of $(J, ip) are given by the equations 
o= ^ cM> = _ 4/2 
dip 

= — = e + 2al + 2bl cos(4<^) + 3c/ 2 (l + cos(4^)) . 



(26) 



As in section [P] , the first of these equations yields (see eq. 

sin(4</?) = , 
COs(4(^) = er = ±1 . 

Thus, the second equation reads 

e + 2(a + ab)l + 3cl 2 (l + a) =0. (27) 

In solving this equation, we shall assume c < 0. 

For a = — 1, the equation is linear. Its only solution reads 

J - = -2(^6) • (28) 

For a = +1, however, we obtain a quadratic equation with two solutions: 



I± = -c~ 1 / 3 [6± Jv + S 2 ) , (29) 



where we have introduced the abbreviations 

^=~6?78' 5 =^73- (30) 

Thus, we obtain three different stationary points corresponding to the three satellite 
orbits in addition to the central periodic orbit at / = 0. 

The dependence of the solutions I on e is shown schematically in figure [| for the 
case \a\ > \b\ and c < 0. Comparing the normal form results to the bifurcation scenario 
described in section |2|, we recognize the sequence of an island-chain-bifurcation at e = 
and a ghost orbit bifurcation at some negative value of e. Thus, our normal form 
correctly describes the bifurcation scenario under consideration, and we adopt it as 
an ansatz function in the uniform approximation (0). This correspondence allows us 
to identify stationary points with classical periodic orbits as follows: The (a = — 1)- 
solution describes the unstable satellite orbits on either side of E c . For E < E' c , the 
(er = +l)-solutions correspond to the asymmetric ghost orbits, whereas for E > E' c , we 
identify the two solutions for a = +1 with the stable satellite orbit (marked by the + 
in figure |5|) and with the additional ghost (— in figure ||). 

We now calculate the stationary values of the normal form, which correspond to 
action differences. They read 



$ ± = 4 (p + 5 2 ) (S ± s[r) + S 2 j + 2 V 5 , 

(31) 

_1 4(o - b) 
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Finally, we need the Hessian determinants of the normal form at the stationary points. 
We calculate them with respect to cartesian coordinates, because in polar coordinates 
the Hessian determinant at the central orbit / = is undefined, and obtain 

Hess ± = [e + 2(o - 3b)I± - 2c/|} [e + 6(a - b)I± + 30c/|} , 

Hess _i = [e + 4aJ_i + Acl^f - tf\ {a - 36 - 2c/_!} 2 , (32) 

Hess o = e 2 . 

We have now found a normal form which can serve as an analytical description of 
the complicated bifurcation scenario. 

5. Determination of normal form parameters 



In order to use the normal form ( |25| ) as an ansatz in the uniform approximation (|TT|), 
we now have to determine the normal form parameters e, a, b, c so that the numerically 
observed action differences are quantitatively reproduced by the normal form results 
(PT|) . Since, according to figure |5], e measures the distance from the period quadrupling 
bifurcation, we choose 

e = E-E c , (33) 



and then solve equations (31) for the parameters a, 6, a 



To achieve this, we introduce 



(34) 



The second equation gives 

r, = h 2 / 3 - 5 2 , (35) 

so that from the first equation we obtain 

5 3 - 3h 2 J 3 5 + 2h + = . (36) 

This is a cubic equation for 5. Its discriminant reads 

D = h\-h 2 _ = y^+§- , (37) 

and, using Cardani's formula, we find its solutions 



^ = \\\- (v^+ v^ 7 ) 2 + y V" (ft*' ^-) 2 ' (38) 

where A G { 1, — \ ± i^ \ is a cube root of unity. If D > 0, A = 1 yields the only real 
solution, whereas for D < all solutions are real. To proceed, we have to choose one of 
the three solutions. 

As can be seen from figure § using the correspondence between stationary points 
and periodic orbits discussed above, we have $ + > 0, and there is an e < such that 
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$_ > for s > e and $_ < for s < e . Consequently, D > for e > e and D < 
for e < Sq. To make 5 real, we therefore have to choose A = 1 if e < £q. 

To determine A for e > Eq, we first observe that the parameters must depend on 
e continuously. Thus, A can only change at energies where equation (|36| ) has a double 
root, viz. D = or e G {0,£o}- Therefore, it suffices to determine A in a neighbourhood 
of e = 0. 

From the plot of action differences in figure |] we find 

$, = a 2 £ 2 + O (^ 



$_ = -r - 0e + O (e 3 ) 



with positive constants a,ft,T. With the help of equations (|38[) and ( pq) we can now 
expand rj in a Taylor series to first order in e: 

ry=fi-(ReA) 2 V 2/3 



4 

2/3 _ XT 4a v^ . 




- (Re A) 2 ) ^ - Re Aim Assign s) £ + O (e 2 ) . (39) 

Requiring this result to reproduce the definition 

1 1 

we find the conditions 

1 , 2aVT 

ReA = --, 3rV3 Slgn£ > ° ' 

which lead to the correct choices of A: 

f 1 : $_ > 

A = 1 .v/3. , • ( 4 °) 

h i signe : $_ < 

I 2 2 5 



Using this result, we can determine r\ and 5 from ( p5[) and (|38[) . Equations (p0|) and 
(Bll) then finally yield the desired parameter values 



c = — 




b = 3c 2 / 3 5 + 



(41) 



8$_i ' 8$_! 

Note that from ( |4"TD the parameters a, 6, and c are explicitly determined as functions of 
the energy e and the three action differences $ + , $_, and $_i. 

In our case, we will determine the normal form coefficients from the scaled action 
differences shown in figure |j. To obtain the actual non-scaled coefficients for different 
values of the magnetic field strength 7, we need to derive scaling laws for the coefficients. 
As we shall display the semiclassical spectra as functions of scaled energy, we prefer not 
to scale the energy difference e = E — E c . Then, with the help of equations (|3lD and 
(J32|), we can convince ourselves that the scaling prescriptions 

a = 7 " 1 / 3 a, 6 = 7 " 1 /3 &) g = 7 -2/3 c ( 42 ) 
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fulfil the requirements of scaling actions according to S — 7 1//3 S' while not scaling Hessian 
determinants. 

6. Evaluation of the uniform approximation 

Now that the ansatz function $ has been completely specified, a suitable approximation 
to the coefficient X in ( |TTD remains to be found. We assume X to be independent 
of <p, and as the value of X is known at the stationary points of / at four different 
values of I (including / = 0), we approximate X by the third order polynomial 
p(I) interpolating between the four given points. This choice ensures that our 
approximation reproduces Gutzwiller's isolated-orbits formula if, sufficiently far away 
from the bifurcations, we evaluate the integral in stationary-phase-approximation. Thus, 
the uniform approximation takes its final form 

d{E) = ^J-^R e exp{j;S (E) - fy} J dY dfy p(I) exp {^>(Y, P Y )} , (43) 

which contains known functions only and can be evaluated numerically. 

It is important to point out, at this juncture, the progress encorporated in equation 
( pE3]) as compared to previous literature results. Starting from the lowest-order normal 
form (PU|), Sieber and Schomerus |7| derived a similar formula for the contribution of 
an isolated generic period-quadrupling bifurcation to the density of states. However, 
as we have seen in section |4.1|, their normal form (|20|) describes the central periodic 
orbit and the stable and unstable satellites only, from which it is evident that the 
uniform approximation given by Sieber and Schomerus cannot take the presence of the 
additional ghost orbit and the occurence of the ghost orbit bifurcation into account. 
The effect of the latter is negligible far above the bifurcation energy E c , where the 
asymptotic behaviour is determined by the three real orbits common to both (Sieber 
and Schomerus's and our own) forms of the uniform approximation. The Sieber and 
Schomerus result, however, is not capable of describing the correct asymptotic behaviour 
below E c , because one of the satellite orbits runs into a bifurcation unforeseen by the 
normal form ( p0|) at E' c , thus causing Sieber and Schomerus's uniform approximation to 
diverge. Below E' c , when the stable satellite orbit does not exist any more, their solution 
also ceases to exist because the required input data is no longer available. Thus, to 
obtain a smooth interpolation between the asymptotic Gutzwiller behaviour on either 
side of the bifurcations, we must make use of the extended uniform approximation fl43|), 
which takes the contribution of the ghost orbit bifurcation into account, as long as we 
do not go to the extreme semiclassical domain where Planck's constant has become so 
small that below E c the asymptotic regime is reached before the ghost orbit bifurcation 
can produce a palpable effect. It is only in this limit that the impact of the ghost orbit 
bifurcation on the semiclassical spectrum vanishes. 

We calculated the uniform approximation (|43"D for three different values of the 
magnetic field strength 7. The results are shown in figure ^. To suppress the highly 
oscillatory contributions originating from the factor exp I ^Sq(E) \, we plot the absolute 
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value of ( fJ5|) instead of the real part. As can be seen, the uniform approximation 
proposed is finite at the bifurcation energies, and, as the distance from the bifurcations 
increases, asymptotically goes over into the results of Gutzwiller's trace formula. Even 
the complicated oscillatory structures in the density of states caused by interferences 
between the contributions from the different real orbits involved at E > E c = —0.342025 
are perfectly reproduced by our uniform approximation (see figures ^|b and []c). We also 
see that the higher the magnetic field strength, the farther away from the bifurcation 
the asymptotic (Gutzwiller) behaviour is acquired. In fact, for the largest field strength 
in figure ^a (7 = 10 -10 ) the asymptotic region is not reached at all in the energy domain 
shown. The magnetic field dependence of the transition into the asymptotic regime can 
be traced back to the fact that, due to the scaling properties of our system, the scaling 
parameter 7 1 ' 3 plays the role of an effective Planck's constant, therefore the lower 7 
becomes, the more accurate the semiclassical approximation will be. 

7. Conclusion 

We have shown that in Hamiltonian systems with mixed regular-chaotic dynamics ghost 
orbit bifurcations can occur besides the bifurcations of real orbits. These are of special 
importance when they appear in the vicinity of bifurcations of real orbits, since they 
turn out to produce signatures in the semiclassical spectra much the same as those 
of the real orbits. Consequently, the traditional theory of uniform approximations for 
bifurcations of real orbits must be extended to also include the effects of bifurcating 
ghost orbits. 

We have illustrated the phenomenon of bifurcating ghost orbits in the 
neighbourhood of bifurcations of real orbits by way of example for the period- 
quadrupling of the balloon orbit in the diamagnetic Kepler problem, and have 
demonstrated how normal form theory can be extended for this case so as to allow 
the unified description of both real and complex bifurcations. 

We picked the example mainly for its simplicity, since (a) the real orbit considered 
is one of the shortest fundamental periodic orbits in the diamagnetic Kepler problem 
and (b) the period-quadrupling is the lowest period-m-tupling possible (m = 4) that 
exhibits the island-chain-bifurcation typical of all higher m. Thus we expect ghost orbit 
bifurcations to appear also for longer-period orbits, and, in particular, in the vicinity of 
all higher period-m-tupling bifurcations of real orbits. 

In fact, a general discussion of the bifurcation scenarios described by the normal 



form (25|) and its more general variant given in |22j for different values of the parameters 
leads us to the conclusion that the appearance of ghost bifurcations in the vicinity 
of bifurcating real orbits is the rule, rather than the exception, in general systems 
with mixed regular-chaotic systems, and thus one of their generic features. It will 
be interesting and rewarding to study higher period-m-tuplings with respect to the 
appearance of ghost orbit bifurcations, and to extend ordinary normal form theory to 
also include the contributions of ghost orbit bifurcations for all higher m. 
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Figure 2. Real and imaginary parts of complex ghost orbits at scaled energy 
E = —0.343. Solid and dotted lines: Stable and unstable ghost satellites created 
at the period-quadrupling bifurcation of the balloon orbit at E c = —0.342025. Dashed 
line: Additional ghost orbit created at the ghost bifurcation at E' c = —0.343605. 
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Figure 3. Ghost orbits at scaled energy E = —0.344 < E' c . Solid and dashed 
lines: Asymmetric ghosts (real parts coincide) created at the ghost bifurcation at 
E' c = —0.343605. Dotted line: Unstable ghost satellite created at the period- 
quadrupling bifurcation of the balloon orbit at E c = —0.342025. 
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Figure 4. Actions, periods and traces of the orbits involved in the bifurcations as 
functions of the scaled energy E = £>y~ 2 / 3 . Solid lines: real orbits, dashed lines: 
ghost orbits symmetric with respect to complex conjugation, dotted lines: asymmetric 
ghost orbits. 
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Figure 5. Sketch of the bifurcation scenario given by the normal form (|25|) for the 
case \a\ > \b\ and c < 0. Solid lines: Real orbits. Dashed lines: Ghost orbits symmetric 
with respect to complex conjugation. Dotted line: A pair of complex conjugate ghosts. 
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Figure 6. Uniform approximation to the contribution of the considered bifurcations 
to the density of states for three different values of the magnetic field strength: (a) 
7 = 10~ 10 , (b) 7 = 1CP 12 , (c) 7 = 10~ 14 . Solid lines: Uniform approximations. 
Dashed lines: Gutzwillcr's trace formula. 



